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ABSTRACT 



The objective of this study was to compare two commonly used 
computer simulation languages: GPSS and FORTRAN. The comparison 
was made by simulating an identical system in both languages. The 
comparison criteria used to evaluate the languages were as follows: 
ability to represent system, simulation time, ability to represent 
stocastic phenomena, programming time, computer running time, monitor- 
ing and debugging, storage requirements, data initialization and starting 
conditions, replication, data collection and display. The results from 
this study indicated that a system may be modeled four to five times 
faster using GPSS as compared to FORTRAN. The modeled system was 
simpler to conceptualize in GPSS, and the model required reduced pro- 
gramming, debugging and monitoring time compared to the FORTRAN 



model . 
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I. INTRODUCTION 



Computer simulation has become an ordinary and useful technique 
of the operations analyst. In industrial and military operations research, 
complex as well as simple systems are being simulated by computer for 
both educational and analytical purposes . In an effort to reduce both 
the amount of detailed programming required and the amount of intricate 
simulation methodology that must be known in order to create a computer 
simulation, several special purpose simulation languages have been made 
available to the operations research community [Refs. 15, 16, 17]. Prior 
to the introduction of these simulation languages, computer simulation 
programs were restricted to either the available assembly or compiler 
languages. It is not clear, however, that the task of simulating a system 
is always simpler when using a simulation language. 

• It is obvious that as a computer language becomes specialized, it also 

becomes restrictive. It is also possible that an increase in speciali- 
zation for a computer language will tend to decrease the efficiency of 
the language in terms of running time and allocated memory space. Moti- 
vated by these thoughts, the objective of this thesis is to compare the 
usefulness of two very common programming languages: GPSS and 
FORTRAN. GPSS is a specialized simulation language for the simulation 
of queueing problems and FORTRAN is a general purpose language. A 
detailed description of these languages is not contained in this thesis 
but may be found in Refs . 3 , 4 , 6 , 7 , 8 , and 14 . 
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The comparison of GPSS and FORTRAN in this thesis is made by 



simulating an identical queueing system in each of these languages. 

For the queueing system simulated, identical assumptions were made 
relative to modeling the system and the same exogenous data was used. 
Models in each of the languages provided comparable output data . At 
the beginning of this study, the author was equally familair with each 
language. Based upon these two simulations, the languages were then 
relatively evaluated according to the following criteria: 

1. Ability to represent system 

2 Simulation time 

3. Ability to represent stocastic phenomena 
4 . Programming time 

5. Computer running time 

6. Data initialization and starting conditions 

7. Storage requirement 

8. Replication 

9. Monitoring and debugging 

10. Data collection and display 

These criteria are discussed in Chapter II in light of their importance 
to simulation model building. 

It should be made clear that the content of this thesis is directly 
related to a comparison of GPSS and FORTRAN in programming a simple 
queueing situation. No attempt is made to compare these languages 
beyond their application in programming this sample problem. 
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II. LANGUAGE COMPARISON CRITERIA 



The criteria used to compare FORTRAN and GPSS represent both the 
economic considerations of choosing a language as well as the desirable 
features of a simulation language. The economic considerations [Ref. 12, 
pgs. 240-241] relate to cost and availability. Cost can be equated to 
the man hours required to obtain a validated computer model, together 
with the computer operating time required to obtain the final results. 
Availability includes the availability of computer hardware as well as the 
availability of knowledgeable programmers and technicians. The desirable 
features of a simulation language [Ref. 9, pgs. 26-38] are qualities that 
assist the programmer in conceptualizing, programming, debugging and 
experimenting with the computer model. A detailed explanation of the 
specific criteria utilized to compare FORTRAN and GPSS follows . 

A. ABILITY TO REPRESENT SYSTEMS 

Most systems of interacting entities can be said to possess two 
structures: static and dynamic. The static structure is independent of 
time, and consists of the framework or structure within which action 
takes place. The dynamic structure is time dependent and includes the 
processes and actions that occur within the static structure of the model. 

In order to simulate a system, it is, therefore, important that the 
language used must contain an inherent framework or flexibility so that 
both the dynamic and static structure of the system can be easily and 
realistically included in the simulation. As examples of the two structures, 
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consider the operation of a message center. The transmitters, receivers, 



and electrical circuits represent the static structure. The arrival and 
departure of messages over the circuits make up the dynamic structure 
for this system . 

B. SIMULATION TIME 

In order to simulate the dynamics of a system, a method is needed 
which portrays the passage of time. Time in a simulation is controlled by 
means of a computer clock which is a programming mechanism that approxi- 
mates the continuous flow of time in the simulation. Clocks are of two 
basic types: fixed time and next event. Fixed time clocks operate by 
advancing simulated time in discrete time intervals. At the completion 
of each time interval, a control mechanism scans all the possible actions 
that may take place and all computation scheduled for this interval is 
performed. The clock then advances into the next interval and the 
scanning process is repeated. Next event clocks utilize a control device 
that, at the completion of an event or computation, advances simulated 
time directly to the time of the next event scheduled to occur. For a 
detailed description of these clock mechanisms see Naylor [Ref. 12]. 

C. ABILITY TO REPRESENT STOCASTIC PHENOMENA 

Events in the real world usually do not occur with regularity or 
according to some fixed pattern. The variability and uncertainty with 
which events occur are most often described by means of probability 
statements and distributions, e.g. , event A has fifty percent chance of 
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occurring, or event B will occur every ten days plus or minus three days. 



To include these probability statements and distributions, a method is 
needed within the dynamic structure of the model to generate random 
numbers. Random numbers used in conjunction with probability distri- 
butions and statements are used to determine the outcome of stocastic 
events, e.g. , did event A occur or not? In order to adequately model 
real world systems, it is therefore necessary to include in the evaluation 
of any computer language the ease with which the variability and un- 
certainty of events can be included in the simulation. 

D. PROGRAMMING TIME AND COMPUTER RUNNING TIME 

For any simulation, programming time and computer running time 
are important considerations. Both times represent dollars and, there- 
fore, must be considered prior to choosing a computer language. Complex 
languages which compile and execute slowly in comparison to simpler 
languages may greatly reduce programming time due to their syntax. In 
deciding on a language to use, a cost and time trade off analysis should 
be made. This analysis can quite often only be made empirically, i.e. , 
by comparing the results of the different languages as is being done in 
this thesis . 

E. MONITORING AND DEBUGGING 

An essential requirement for any programming language is that it 
ought to provide features that assist the programmer in debugging syntax 
and logic errors. For simulation languages in general, three desirable 
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debugging features are: 



1. A list of compilation and execution errors referenced by 
statement number and clock time. 

■ 2 . A list of execution errors reported with a complete printout 

of the status of the program at the time of error. This print- 
out should contain as a minimum: clock time, status of 
appropriate variables and related parameters from the sub- 
routines currently in effect. 

3 . A trace back feature which would provide some method to 
trace a given event thru the system. 

These three features are also necessary in order to monitor the system 
dynamics. Without the ability to trace an event and the ability to 
observe the system at various times, many logic errors would be im- 
possible to locate [Ref. 10, pgs . 35-36]. 

F. DATA INITIALIZATION AND STARTING CONDITIONS 

Simulation studies often use extensive exogenous data. Therefore, 
convenient methods to initialize exogenous data into the computer model 
is a desirable language feature. Once the model has been initialized, 
the immediate objective then becomes to obtain information about the 
system at various points in time. In the analysis of some systems, the 
analyst may decide to study the system during steady state conditions, 
while for other systems, the transient conditions may be more important. 
Another desirable feature for comparison of language is then the ease 
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with which the inherent structure of the language allows the analyst to 
drive the simulation to a steady state condition. This must include, of 
course, a convenient method for minimizing the amount of execution time 
required to arrive at a steady state condition. 

G. STORAGE REQUIREMENTS 

Computer core storage is fixed for each model computer. This re- 
stricts the capability of each model computer to perform jobs within its 
memory capacity. In order to obtain the most efficient use of a computer's 
memory capacity, many computer facilities encourage programs to be 
written with small storage and time requirements by providing quicker 
turn around times to these programs. The language storage requirement 
therefore becomes important in obtaining shorter turn around times during 
the debugging and experimentation phase of the simulation. Core storage 
requirements of a computer language can therefore materially affect the 
time required to project completion. 

H. REPLICATION 

When simulation results are put to the test of statistical analysis, 
it is of primary interest to obtain results with a specified degree of 
confidence. In order to be able to specify a degree of confidence, a 
number of observations of the experiment must be obtained. These 
observations are obtained by replicating the computer experiment the 
desired number of times. The ease with which program replications can 
be made is therefore a measure of contrasting two simulation languages. 
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I. DATA COLLECTION AND DISPLAY 

Data collection is a necessary requirement for any simulation. 
Ideally, data collection should be automatic, but in any event, data 
collection statements should not interfere with the operations of the model. 
Since completely automatic data collection is not practical, the next best 
alternative is to have certain information automatically collected and the 
remaining information available at the discretion of the programmer. The 
following data should be available as output at little or no inconvenience 
to the programmer [Ref. 10, pgs . 32-34]. 

1 . The number of observations and the maxima and minima for 
all variables 

2. Sums and sums of squares for time independent variables 

3 . Time-weighted sums and sums of squares for time-independent 
variables 

4. Variable value histograms for time-independent variables 

5. Time-in-state histograms for time-dependent variables 

6. Time series plots over specified time intervals 

The recording of results is a primary objective in every study. 
Therefore, the language selected for simulation should allow for varying 
output formats that are dependent upon the programmers needs. However, 
as is the case with data collection statements, the programmer should not 
have to spend a great deal of time in writing or debugging output data 
statements . 
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III. THE MODEL 



This chapter introduces the torn tape relay network concept of 
operation in Section A and describes the system that was simulated for 
this case study in Section B. The FORTRAN and GPSS models of the torn 
tape relay network are described in Sections C and D respectively. The 
FORTRAN and GPSS programs are discussed with the intent of providing 
the reader insight into the methods that can be used by a programmer to 
simulate the static and dynamic structure of a system. 

A. TORN TAPE RELAY CONCEPT 

The purpose of a torn tape relay network is to move message traffic 
from an origin to a destination within the network. A typical torn tape 
relay system consists of a number of terminal and relay stations con- 
nected by circuits. A terminal originates and terminates messages for a 
military headquarters. Messages originated at a terminal are translated 
into a coded perforated tape. The perforated tape code is then transmitted 
to a relay station according to the assigned precedence and desired routing 
of the message. A relay station is typically connected to many terminal 
stations in the geographical area as well as to other relay stations. The 
mission of a relay station is to accept messages from the supported 
terminals and connecting relays and to retransmit the messages to the 
next station according to the messages routine instructions. This next 
station may be another relay or a supported terminal station. 
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A message may be assigned one of four precedences. Listed in 
order of importance from high to low, these precedences are: flash, 
immediate, priority, and routine. All messages are transmitted from one 
station to another on a first in first out, FIFO, discipline ordered by 
precedence. When a flash message is received at a station and no 
circuit is free to the next station in accordance with the routing instruct- 
ions, a lower precedence message occupying a transmitting circuit is 
p re _ em pted. The flash message is then transmitted over the pre-empted 
circuit and the pre-empted message is refiled according to its original 
receipt time. The pre-empted message is then eventually retransmitted 
according to the normal FIFO discipline. A message with a flash pre- 
cedence is the only type of message with pre-empting authority. 

Prior to reaching the final destination, messages often remain in 
queue at some station and are said to be backlogged. Station backlogs 
are computed to be the time required to transmit all messages in queue 
for a given station. The total backlog is often separated into two 
categories: a high precedence message backlog consisting of flash and 

immediate messages and a low precedence backlog consisting of routine 
and priority messages. Backlog statistics are used as a measure of 
effectiveness to evaluate a station's ability to pass traffic. 

The number of circuits interconnecting two stations regulates the 
backlog that exists between these stations. Additional tape transmitters, 
receivers, and personnel to operate the equipment are also needed as the 
number of circuits increase. The additional circuits, equipment and 
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personnel can be equated to dollars. Hence, in designing or operating a 



torn tape relay, one attempts to minimize the cost of operation by mini- 
mizing the number of interconnecting circuits, subject to some fixed 
effectiveness criterion. The maximum backlog that will be tolerated is 
often used as the fixed effectiveness criterion. 

B. SYSTEM MODELED 

The torn tape relay network that was modeled for this study con- 
sisted of five terminal stations and one relay station. A diagram of the 
system simulated is depicted in Figure 1. The terminal stations are 
numbered 1, 2, 3, 4 and 5. The relay consisted of five bays and are 
labeled 6, 7, 8, 9 and 10 in the diagram. Each bay contained relay 
equipment used to transmit and receive messages to a given terminal. 

The circuits interconnecting the terminals and relay bays are represented 
by directed lines in the diagram. Transmitting circuits are indicated by 
the direction of the arrow. 

The five terminal stations and five relay transmit bays together 
with the associated circuitry comprised the static structure of the system. 
The dynamic structure of the system consists of messages possessing a 
given precedence, arrival time, message length and destination moving 
thru the static structure. 

The problem addressed in this study was to simulate the torn tape 
relay network represented in Figure 1 and to determine the minimum 
number of circuits required between stations consistent with acceptable 
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TORN TAPE RELAY NETWORK 




Figure 1 



backlog figures. Backlogs were considered excessive if low precedence 



messages were backlogged in excess of 15 minutes and high precedence 
messages in excess of 2 minutes . 

Various operating characteristics for all the stations were assumed 
for the system. These characteristics included probability distributions, 
message density functions, peak load figures and average message 
length. The following assumptions were used in modeling the torn tape 
relay network: 

1. All messages introduced into the system reached the desired 
destination and no messages were misrouted. 

2 . Message center record keeping procedures did not contribute 

to backlog figures . 

3 Circuit outages were not considered; hence, all circuits 
experienced one hundred percent reliability. 

4 . The time between message arrivals for each precedence 

message was distributed exponentially with mean peak hour 
values as listed in Figure 11, Appendix A. 

5. The peak hour arrival rate was dependent upon time of day 
and was different for each site. The peak hour arrival rates 
utilized are listed in Figures 5 thru 9, Appendix A. 

6. The average message length was not time dependent, but was 
dependent upon the station and precedence of the message. 
Message length was assumed exponential with mean values 
as listed in Figure 12, Appendix A. 
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C. FORTRAN SYSTEM MODEL 

FORTRAN is a general purpose language, GPL, in that the language 
was not designed for special purpose applications such as simulation. 
FORTRAN is a comparatively simple language to learn, because the symbols 
and mnemonics are analogous to those of mathematics. For this reason, 
and due to the FORTRAN compilers universal availability, the language 
became a common computer language for scientific computations. 

FORTRAN makes use of variables, subscripted variables, constants, 
expressions and functions. The number of different functions or sub- 
routines used by a programmer is limited only by the size of the computer 
available. For simulation programming using FORTRAN, the programmer 
has great flexibility by being able to write a general main program con- 
taining a timing clock, and then calling various subroutines for execution. 

In order to model the static structure of a system in FORTRAN, one 
must portray the structure of the system using the subscripted-unsub- 
scripted variables and constants in conjunction with FORTRAN functions 
and expressions. The dynamic structure of a system is represented in 
FORTRAN by altering the flow of statement execution in a manner that 
reflects the action occurring in the actual system. The order of state- 
ment execution is varied in FORTRAN by use of IF statements, COMPUTED 
GO TO statements and logical type statements. 

In order to simulate the torn tape relay network in FORTRAN, the 
static structure was represented by subscripted arrays. Four three 
dimensional arrays, namely RR (i,j,k), PP (i,j,k), OO (i,j,k) and 
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FLASH PRECEDENCE ARRAY FF(i,j,k) 
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Figure 2 



FF (i,j,k), were used to store information pertaining to routine, priority, 
immediate and flash messages, respectively. Figure 2 contains a repre- 
sentation of the FF array. The k dimension of each array was equal to 
three, and was used to store the arrival time, required transmission 
time, and destination of each message waiting for transmission at the 
terminals and relay transmit bays. The i and j dimension of these arrays 
fixed respectively the number of messages and the station at which the 
messages were currently located. The j dimension was equal to ten with 
argument one thru five corresponding to the five terminals and arguments 
six thru ten corresponding to the five relay transmit bays. The arguments 
of the j dimension corresponded to the station numbers depicted in 
Figure 1 . 

The d eimension was different for each precedence array and was 
selected to permit at least a thirty minute backlog of messages to exist 
for each precedence category at each station. For k equal one, message 
arrival times were stored in the i dimension in increasing order. 

Ten one dimensional arrays HOLDl(i) thru HOLDlO(i) were used to 
represent the circuits between terminals and relay bays. The arrays 
HOLD1 thru HOLDIO correspond to the numbering of Figure 1; these 
arrays contained the time at which a specified circuit would be free. 

The number of transmit circuits from one location to another was equal to 
the dimension of the HOLD array corresponding to this pair of stations. 

The dynamic structure of the system was represented by moving 
message parameters from the terminal storage areas represented by 
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position 1 thru 5 of the j dimension to the desired relay transmit bay 
represented by position 6 thru 10 and then from the relay to the desired 
destination. The order of statement execution was controlled by a FORTRAN 
COMPUTED GO TO statement that acted as a next event clock. Program 
control was transferred to the station that had the highest precedence 
message with the lowest arrival time available for transmission. A 
message was available for transmission if the scheduled message arrival 
time was less than or equal to the simulation clock time and there was a 
free circuit in order to transmit the message to the connecting station. 
Transmission from a terminal was accomplished by transferring message 
parameters from the array column representing the terminal to the array 
column representing the relay transmit bay that will transmit the message 
to its final destination. In addition, the message transmission time is 
added to current simulation time and this value is placed in the proper 
HOLD array. 

For example, suppose a flash message from terminal two to terminal 
three was selected as the next event in the simulation. Since relay 
transmit bay eight transmits to terminal three, the message parameters 
from column two would be transferred to column 8 of array FF. The 
arrival time of the message for terminal eight would also be updated by 
the transmission time of the message. In addition, the time at which 
transmission ends would be transferred to HOLD2(i). 

The information was deleted from the HOLD array after the required 
message lapse time by use of the next event clock. Each time a set of 
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message parameters was transferred from a terminal column, the remaining 



message parameters moved forward one position, and a new set of message 
parameters was created to occupy the empty position. After the creation 
of a new message, all positions in the j dimension column were filled, 
and the message arrival times were in ascending order. Once a message 
was transferred to a relay transmit bay column, the message arrival 
times in the column were not necessarily in ascending order. Hence, 
after a message was transferred to a relay transmit bay column, all 
message parameters in the column were ordered to insure the lowest 
arrival time message was in position one and therefore was the next 
scheduled message to be transferred from the column. 

Messages in a relay transmit bay column selected as the next event 
were processed in a similar manner. The single exception being that 
messages transmitted by the relay were destined for the final destination. 
Hence, after a message was held in the respective HOLD array for the 
messages transmit time, the message attributes were destroyed. The 
destruction of the message's parameters represented the successful 
delivery of the message. 

Backlogs were computed after transferring a message to the next 
station. A message was backlogged if simulation time was greater than 
the message arrival time stored in the j dimension of the respective 
precedence array. All messages stored in relay terminal bay columns 
were backlogged, whereas, only the messages with arrival times less 
than clock time were backlogged in the terminal columns. The messages 
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with times greater than the clock time in the relay columns were future 
events and had not yet been introduced to the model. The backlog for a 
station was computed by totaling message time lengths for those messages 
with arrival time less than or equal to clock time. The maximum backlog 
that occurred during a computer run was recorded for each terminal station 
and relay transmit bay. 

A flow chart of the torn tape relay network in FORTRAN logic is 
illustrated in Figure 3. This is an abbreviated flow chart and includes 
only the main program and the programmed operation of a terminal. Al- 
though some detail is omitted, Figure 3 portrays a system flow chart 
for FORTRAN. A complete flow chart of the FORTRAN model is included 
in Appendix B, together with an explanation of the symbols, and descrip- 
tion of the arrays and subprograms used. 

D. GPSS SYSTEM MODEL 

General Purpose Simulation System, GPSS, is a popular simulation 
programming language. GPSS was designed to be used as a simulation 
language and as such makes special features available to simulation 
programmers . These features hopefully reduce the time required to 
simulate a system by reducing problem formulation time, programming 
time and debugging time. GPSS was designed by the International 
Business Machines Corporation and is well documented. 

GPSS uses a block diagram flow chart procedure to represent the 
static structure of a system. The blocks used to represent the static 
structure are divided into six categories. The most fundamental of 



30 



these categories contain the GPSS entities STORAGES, FACILITIES and 



LOGIC SWITCHES. These entities are the building blocks of GPSS. A 
FACILITY allows entry to only one transaction at a time. A transaction is 
a dynamic entity and is used to represent movement in the system. The 
length of time a FACILITY is occupied by a transaction may be determined 
by a transaction's parameter value. A parameter describes a particular 
transaction attribute. STORAGES allow entry to several transactions 
simultaneously. LOGIC SWITCHES are two position gates which alter 
transaction flow according to the status of some other entity or trans- 
action parameter value. The dynamic structure of a system is modeled 
in GPSS by the use of transactions . Transactions are created and 
destroyed as needed during the computer simulation. Transactions may 
have associated with them a set of parameters. The values assigned to 
these parameters remain with the transaction until they are changed or 
the transaction is deleted from the model. 

Associated with GPSS entities are certain standard numerical 
attributes. The value of these attributes may be referred to or collected 
by the programmer during the execution of the program . Standard 
numerical attribute values may be used to trigger LOGIC SWITCHES, 
assigned to transaction parameters, as well as stored for future use. 

For the torn tape relay network, the terminal and relay sites, 
transmitters, receivers, and interconnecting circuits were modeled by 
use of a block diagram. The messages were represented by transactions 
which moved through the block diagram according to the instructions of 
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the model. Transactions were created in the model for each precedence, 
and from each terminal according to an exponential distribution. The 
intermediate and final destinations, precedence, and message length 
were assigned to transaction parameters in the torn tape relay model. 

Extensive use was made of the FACILITY entity. FACILITIES were 
used to represent torn tape transmitters. Transactions waited in and 
advanced in queue according to the FIFO discipline. A transaction 
occupied a facility, for the simulation time length of the message. A 
PREEMPT block allowed transactions representing flash messages to 
pre empt a FACILITY occupied by a lower precedence transactions. 

Entities may be referred to by indirect addressing in the GPSS 
program. Normally, each entity is designated by number or mnemonic 
in the flow chart of the system. However, it is also possible to 
designate a particular entity by the value of a transaction parameter. 
Indirect addressing is often used to specify entity numbers when there 
are several identical processes occurring in the system to be modeled. 

In the torn tape relay model, five terminals and five relay trans- 
mitter bays performed identical operations. Each location transmitted 
messages according to the FIFO discipline. The terminals transmitted 
to the relay, and each relay transmitter bay transmitted messages to a 
particular terminal. Therefore, since the functions at each location 
were identical, the required programming was reduced tenfold by in- 
corporating indirect routing for referencing most entity blocks. 
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A flow chart in GPSS of a single terminal is depicted in Figure 4. 



This flow chart differs slightly from the actual GPSS model used in that 
indirect routing, statistic gathering and storage of transactions in user 
chains is omitted. A complete flow chart of the model is included in 
Appendix C together with an explanation of symbols used. A comprehensive 
discussion on GPSS may be found in IBM manuals and Schriber [Refs. 6-8,14]. 
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IV. GPSS-FORTRAN COMPARISON 



GPSS and FORTRAN are evaluated in this chapter using the criteria 
discussed in Chapter II. The evaluation is based upon observations and 
experiences encountered while modeling the torn tape relay network dis- 
cussed in Chapter III. 

A. ABILITY TO REPRESENT SYSTEM 
1 . FORTRAN 

The torn tape relay network was first simulated in GPSS and 
then in FORTRAN. The author was hence quite familair with the system 
prior to programming the system in FORTRAN. Even so, a considerable 
amount of time was consumed trying to conceptualize the problem in the 
form of subscripted-unsubscripted variables . Many methods proved 
infeasible due to the requirement to record backlog statistics and due 
to the pre-empt authority of flash messages. 

In order to collect backlog statictics , backlogged messages 
were simulated by creating messages in advance with a scheduled time 
of arrival. When simulation time advanced and equaled the arrival time 
of a message, if a circuit was available to the desired destination then 
the message was transmitted. If a circuit was not free, the message 
would wait in queue and be processed according to the FIFO discipline. 
The only practical method to store messages waiting for an available 
circuit, appeared to be to store the message parameters in arrays. Since 
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arrays were used, the following decisions had to be made before a system 



static structure could be formulated: 

1. How should message parameter values be assigned to 
arrays ? 

2. How many arrays and what dimensionality should be 
used ? 

3. What properties for each message should be retained? 

4. How is information to be exchanged between arrays 
in order to represent the dynamic structure of the 
system? 

2 . GPSS 

The torn tape relay network was comparatively simple to 
represent in GPSS. The system static structure was logically represented 
by means of a GPSS block diagram similar to Figure 4 of Chapter III. The 
blocks FACILITY, STORAGE, TRANSFER, QUEUE and PREEMPT readily 
portrayed the realistic operation of the system. Once the entire block 
diagram was completed, the programming of the model in GPSS followed 
directly, since there exists a one to one correspondence between GPSS 
blocks and GPSS program statements. The dynamic structure of the torn 
tape relay problem was represented in GPSS by GPSS transactions. For 
this problem, transactions were considered messages. Transactions 
were created at prescribed interarrival rates and assigned parameter 
values corresponding to precedence, message length and destination. 

The flow of transactions followed the block diagram structure of the 
network . 
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3 . 



Comparison 



The torn tape relay network system was much easier to 
represent in GPSS than in FORTRAN. The construction of the GPSS block 
diagram using the block mnemonics facilitated building a GPSS model of 
the system. Initially the system was difficult to conceptualize using 
FORTRAN statements since it was possible to become involved in the 
details of FORTRAN and lose presence of the system to be modeled. The 
flexibility of modeling a system in FORTRAN can be a handicap for less 
experienced programmers because a programmer must select from among 
apparent alternate methods of structuring the model. Hence, it can be 
difficult and time consuming to determine which methods are feasible 
and which feasible method is best. 

B. SIMULATION TIME 
1 . FORTRAN 

The incremental time clock and the next event clock were 
considered for use in the FORTRAN model. The incremental time clock 
was not used, however; this clock appeared to be simple to utilize but 
did not appear very efficient since simulation time would be advanced 
in increments of one second and many hours of operation had to be 
simulated. The next event clock was used and appeared to be more 
efficient since simulation time is advanced directly to the time of the 
next event . 

In retrospect, the Incremental clock might have been a more ef- 
ficient choice, since time progressed in small time increments in the 
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model. If an incremental clock had been used, FORTRAN bolean variables 
could have been incorporated in lieu of the FORTRAN IF statements and 
the COMPUTED GO TO statement. Bolean variables associated with each 
possible event could have been set to a go condition if the event was 
scheduled to occur during the next time increment. Bolean variables 
require considerably less computer time and hence might have resulted 
in a more efficient time clock. 

2 . GPSS 

The timing mechanism of GPSS was a next event clock. The 
GPSS program operated by moving transactions thru a block structure 
that represented the static structure of the system. Each block repre- 
sented the possible occurrence of an event in the system. The GPSS 
master program maintained a record of when these events were scheduled 
to occur and processed the events in their proper sequence in time. 

3 . Comparison 

The most important difference in comparing GPSS and FORTRAN 
with respect to simulation time is that GPSS has a built in timing clock 
and FORTRAN does not. When using GPSS the simulation of a system 
must include the next event clock. However, no programming or struct- 
uring of this clock mechanism is required since the clock feature of 
GPSS is entirely automatic within the block structure. The simulation 
of a system in FORTRAN does require the programming and structuring of 
a clock mechanism to meet the specific needs of the system being 
simulated . 
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C. ABILITY TO REPRESENT STOCASTIC PHENOMENA 

1 . FORTRAN 

The mathematical structure of FORTRAN facilitated the pro- 
gramming of probability distributions, functions and functional relation- 
ships. Many useful mathematical functions are available directly from 
the FORTRAN library. In addition, subroutines are available to generate 
random numbers and variates from most of the ordinary theoretical 
probability distributions. Variates can also be easily generated from 
empirical distributions , see Naylor [Ref. 12]. 

Random numbers were generated in the FORTRAN torn tape 
relay model using a library subroutine. The inverse transformation 
technique was then used to obtain variates from an exponential proba- 
bility distribution which required the evaluation of a natural logarithm. 
Several methods were available to generate random numbers and expo- 
nential variates, some of which were more efficient than those used. 

The methods used were selected due to their simplicity and immediate 
availability . 

2 . GPSS 

GPSS contains eight independent random number generators. 
These generators can be synchronized or can operate independently at 
the discretion of the programmer. 

Two methods are available to generate variates from proba- 
bility distributions other than uniform. One method uses the GPSS 

( 

sta tement and a CDF approximated by linea r segments . This 
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method generates variates by using a random number argument and linear 
interpolation on the CDF approximation. This method was used to 
generate exponential variates for the torn tape relay model. 

The second method uses a GPSS VARIABLE statement to express 
the CDF in closed form. Variates can then be generated using the CDF's 
inverse transformation function or some other acceptable method. This 
method is inconvenient, however, since GPSS contains no internal 
routines to evaluate functions such as natural logarithm, sine, gamma, 
squareroot, etc. These functions must be evaluated by using GPSS 
VARIABLE statements to include the evaluation of power series, etc. 

3 . Comparison 

GPSS provides adequate methods to generate variates from a 
probability distribution in order to represent stocastic phenomena . The 
presence of eight random number generators was convenient, but the 
lack of GPSS routines to compute common mathematical functions could 
prove inconvenient. Stocastic phenomena are easily represented in 
FORTRAN and subroutines are available to generate random numbers. 

Any number of independent random number generators can be included in 
the FORTRAN program. Variates from probability distributions can be 
obtained by linear interpolation of the approximated CDF or from the 
theoretical distribution in both languages. However, the availability 
of FORTRAN library functions often simplifies computing closed form 
probability functions compared to the GPSS method. 
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D. PROGRAMMING TIME 



1 . FORTRAN 

FORTRAN statements can be written with or without the aid of 
a flow chart. A general flow chart containing narrative descriptions of 
the actions that are to be portrayed in FORTRAN is often quite useful. 
Such a flow chart helps solidify model concepts and often suggests 
methods of approach. A detailed statement by statement flow chart is 
seldom required except for extremely complicated portions of a program. 
The writing of the FORTRAN program for the torn tape relay problem, 
after a narrative flow chart had been developed, was comparatively 
simple. The FORTRAN program required in excess of 600 statements to 
model the system. 

2. GPSS 

GPSS programming time was considered one of the major 
advantages of the language. The GPSS block diagram used to simulate 
the system greatly simplified the statement programming of the model. 
Once the block diagram had been created, it was a simple matter to 
write the program statement associated with each GPSS block. The 
GPSS program consisted of 200 statements of which fifty were FUNCTION 
and VARIABLE statements and sixty were used to collect statistics for 
output purposes . 

3 . Comparison 

The programming phase of the simulation in GPSS was sig- 
nificantly faster than in FORTRAN. Often where one statement was 
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needed in GPSS to represent an event, a subroutine was required in FORTRAN 
to represent the same event. Consequently, the FORTRAN program required 
three times as many statements as needed by GPSS in order to perform 
the same task. The time required to program the FORTRAN model was 
significantly higher than the three to one ratio indicated by the number 
of statements required. It is quite conceivable that a GPSS programmer 
could be experimenting with the completed GPSS model while a FORTRAN 
programmer was still investigating a method of attack to program the 
same system . 

E. COMPUTER RUNNING TIME 
1 . FORTRAN 

The FORTRAN model required seventeen minutes to compile 
and execute the same number of runs conducted by the GPSS model in 
six minutes. One would generally expect that the FORTRAN program 
would execute faster than the GPSS program. This expectation is based 
on the fact that FORTRAN is a lower level language and utilizes less 
sophisticated data structures and methods to vary the flow of statement 
execution. An attempt was made to determine the reason for this lengthy 
execution time . 

The probable cause of the excessive time was due to the 
method used to advance messages in arrays RR, PP, OO, FF . After a 
message was transferred from a terminal station array column, all 
message attributes were moved forward one position prior to a new 



43 



message being created. After the message was transferred to a relay 
transit bay column, all message arrival times were ordered to insure the 
message with the smallest arrival time was in position one of its respective 
column. The transmission of one routine message required changing 
several hundred computer storage values. This perpetual changing of 
storage values undoubtedly resulted in a high FORTRAN program execution 
time . 

2 . GPSS 

The GPSS model of the torn tape relay network required six 
minutes and five seconds to execute eight repetitions of the torn tape 
networks peak period of operation. 

3 . Comparison 

Although it may be possible to write faster executing programs 
in FORTRAN than in GPSS, the torn tape relay example indicates that 
this is not automatic. It appears that the GPSS language includes 
efficient routines which when assembled by a programmer into a computer 
program provide a comparatively efficient method of building a model. 

A GPSS program may not be as efficient as theoretically possible in 
FORTRAN. However, for the inexperienced programmer, the efficiency 
of the GPSS language may be difficult to match. 

The GPSS program written for this study was three times 
faster in computer execution time than the FORTRAN program. Twenty- 
five replications could have been obtained using the GPSS program as 
compared to eight using the FORTRAN program in the same amount of 
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computer time. These seventeen additional observations would have sig- 



nificantly increased the sample size and provided more confidence in the 
results obtained. 

F . MONITORING AND DEBUGGING 
1 . FORTRAN 

No set debugging scheme was available to ferret out logic 
errors in the FORTRAN program. The syntax diagnostics available are 
quite often very hard to correlate with program logic errors. As an 
example the FORTRAN program developed several undesired loops that 
held simulation time constant while exhausting computer execution 
time. The location of these errors was difficult to determine in a single 
diagnostic run since the simulation clock time of the errors had to be 
determined. This was done by printing the simulation time after each 
message was processed. The printout was then used to determine the 
simulation time when the undesired loop first occurred. A second diag- 
nostic run was then used to obtain information dumps corresponding to 
the model simulation time encompassing the undesired loop. A dump 
consisted of the simulation time and the relevant arrays. By using 
this data, logic errors were isolated by tracing the movement of message 
parameters among array columns. This method was quite time consuming 
since there existed in excess of 3,000 relevant storage locations. On 
several occasions the error could only be isolated to a particular sub- 
routine and additional runs and output were required to locate the error. 
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Once the location of the undesired loop was located, however, it was a 



comparatively simple task to determine the cause and correct the error. 

2 . GPSS 

GPSS debugging and monitoring facilities were excellent. 
Chapter 17 of the IBM User's Manual [Ref. 7] provides the programmer 
with many valuable hints on monitoring and debugging the GPSS program. 
Errors detected during execution of the program triggered an elaborate 
information dump which automatically provided the following information: 
1. A coded error message 

2 . Simulation time error was encountered 

3 . The block number a transaction was currently in 

4. The next block a transaction was to be processed by 

5. The number of transactions processed by each block 
in the program 

6. Current and future event chains 

7. GPSS FACILITY statistics 

8. QUEUE statistics 

9. STORAGE statistics 

This dump included ample information to isolate the majority of the 
errors encountered debugging the torn tape relay model. However, on 
occasion, an additional monitoring and debugging feature was required. 
This feature was the GPSS TRACE option. The TRACE feature provided a 
capability of tracing a transaction with a given attribute thru the static 
structure of the model. Transaction location and parameter values were 
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printed as the transaction followed the GPSS block diagram. In addition, 
information dumps were easily obtained at strategic locations in the pro- 
gram by using a traced transaction to trigger the dump. The dynamics of 
the model were then reconstructed using both the information dumps and 
the trace block data. From this information, the cause of error was 
easily located. 

3 . Comparison 

The GPSS automatic information dump that occurred when an 
error was detected facilitated the isolation of 90 percent of the errors 
encountered. When additional diagnostic information was needed, the 
GPSS TRACE option was available. These two features greatly simplified 
the task of monitoring the debugging the GPSS torn tape relay network 
program. In contrast, the FORTRAN debugging and monitoring error 
detection schemes were not suited for detecting logic errors in the simu- 
lation model. The method used by the FORTRAN programmer must depend 
on the type of error encountered and the particulars of the model. Error 
detection required ingenuity, time and often luck in order to locate the 
cause of errors in the FORTRAN torn tape relay model. 

G. INPUT DATA AND INITIALIZATION 
1. FORTRAN 

Exogenous data for the torn tape relay was made available 
to the model conveniently by use of FORTRAN DATA statements . In 
particular it was found by experimentation that the model did not have 
to be preloaded with backlogged messages at the beginning of the first 
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run. Using DATA statements, the terminal stations and relay transmit 
bay stations corresponding to the columns of arrays RR, PP, OO , FF 
could have been pre-filled with a desired number of backlogged messages 
of any length, destination and arrival time. 

2 . GPSS 

Exogenous data for the GPSS version of the torn tape relay 
network was initialized into the model by means of the GPSS FUNCTION. 
Although transactions were not preloaded into the GPSS version of the 
model to represent a backlogged message condition, preloading could 
have been easily accomplished using the GENERATE, and MATRIX blocks. 
In GPSS a programmer can specify an upper limit on the number of trans- 
actions that will be created by a particular GENERATE block and GENERATE 
blocks could have been used to create the desired number of backlogged 
messages for each precedence category and station. 

3 . Comparison 

Exogenous data initialization was handled adequately by 
both GPSS and FORTRAN for the torn tape relay network system. Although 
the model was not preloaded with backlogged messages for start-up, it 
appeared that both GPSS and FORTRAN could have handled the requirement 
without difficulty. 

Although DATA statements were used for exogenous data in- 
put for the FORTRAN model, additional methods were also available, 
e.g. , data could also have been inputed by means of magnetic tape or 
punched cards. GPSS relies heavily on the FUNCTION statement, 
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SAVE VALUE, MATRIX and INITIAL blocks to handle exogenous data require- 



ments . GPSS cannot accept magnetic tape input and can only read data 
cards with the assistance of the HELP block. However, the IBM manual 
[Ref. 7], states that the HELP block is intended only for users who are 
thoroughly familiar with the internal operation of the GPSS program. This 
last requirement tends to eliminate many hopeful simulators. 

GPSS and FORTRAN adequately handled the exogenous data 
requirement for the torn tape relay network. However, the FORTRAN 
language would be preferable to GPSS if the system to be simulated 
required voluminous input data . 

H. STORAGE REQUIREMENTS 
1 . FORTRAN 

The FORTRAN model of torn tape relay network required 64K 
bytes of computer storage. 

2. GPSS 

The GPSS model of the torn tape relay network required 12 8K 
bytes of computer core storage. GPSS can be operated on the IBM/3 60/67 
in either a 128K or 256K mode. The 256K mode enables a programmer to 
utilize an increased number of entities of each category. The 128K 
version allocates 150 STORAGE blocks, whereas the 256K version 
allocated 300 STORAGE blocks . The only critical factor is the total 
amount of storage required by a program. Unused space from one entity 
group can be reallocated to another entity group and this allows for an 
increase in the number of blocks allocated to a specified entity group. 
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An area referred to as GPSS COMMON also requires allocation of storage 



space. GPSS COMMON is the space used dynamically for transaction 
information, temporary data, and for storing statistical data. The GPSS 
torn tape relay network model required COMMON storage in excess of 
that allocated to the 12 8K version. The additional storage for the torn 
tape relay model was obtained from unused entity blocks. 

3 . Comparison 

The FORTRAN model required half the storage required by 
the GPSS model, that is 64K versus 128K. This reduced storage require- 
ment did not reflect a faster turn around time for the FORTRAN program 
due to the long execution time of the program. 

I . REPLICATION 
1 . FORTRAN 

Replication for the FORTRAN version of the torn tape relay 
network was accomplished by use of a FORTRAN IF statement in con- 
junction with the next event clock. Prior to selecting the next event, 
a test was conducted to determine if the current simulation time was 
in excess of the maximum simulation time for a run (9,000 seconds). 

If simulation time was in excess of 9,000 seconds, statement execution 
was transferred to record pertinent statistics. After the statistics were 
recorded, the number of replications completed was incremented and 
tested. If the required number of runs had been completed, the simu- 
lation was terminated, if not, the clock was reset and another repli- 
cation was performed. 
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GPSS 



2 . 

Replication is accomplished in GPSS by use of RESET and 
CLEAR blocks and by the redefinition of blocks . The RESET block sets 
all statistical tables to zero except those specified by the programmer. 
The RESET block does not alter the status of any block entities or alter 
the values of any transaction parameter values. The CLEAR block per- 
forms a function similar to the RESET clock except that in addition all 
transactions currently in the system are destroyed. 

The GPSS torn tape relay model was replicated with the aid 
of RESET blocks and the redefinition of entity blocks in conjunction 
with a specially constructed timing program. The timing program was 
a simple subprogram that controlled the simulation time for each rep- 
lication. The program consisted of a GENERATE block that created a 
transaction. This transaction then waited in an ADVANCE block for the 
desired time length of each replication (9,000 seconds). The trans- 
action was then used to gather and store desired model statistics in 
GPSS TABLES. After the statistics were gathered, the transaction was 
destroyed. A RESET block then caused entity block statistics to be 
zeroed except for specified tables. A START block was then redefined 
which caused a new transaction to be created and the above process 
was repeated . 

3 . Comparison 

The GPSS features available to assist the programmer in 
replication are very powerful. Specific entity blocks were designed 
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to assist the programmer in every facet of replication. Undoubtedly, one 
could program the same features into a FORTRAN program, but to do so 
would be programmer time consuming and require programmer ingenuity. 

J. DATA COLLECTION AND DISPLAY 

1 . FORTRAN 

FORTRAN does not provide any output automatically and all 
output desired must be programmed for within the model. This may often 
require ingenuity and a considerable amount of programming time. Of 
importance , however , is that desired output can always be obtained in 
any format. 

2 . GPSS 

The GPSS program stored statistics in GPSS SAVEVALUES 
during the execution of a replication, and the GPSS timing program was 
used to collect desired statistics after each replication. The statistics 
were stored until the completion of the last replication and then printed 
in tabular form . 

Some output from GPSS is automatic while other output can 
be programmer specified. The automatic output consists of the following: 

1 . Block counts 

2 . Current value stored in SAVEVALUES 

3 . USER CHAIN statistics 

4 . QUEUE statistics 

5. FACILITY statistics 
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6 . 



STORAGE statistics 



7. Relative and absolute clock times 
The automatic output is voluminous and provides general information about 
the system . 

Information relating to transaction,parameters time and specific 
facets of the system must be programmer specified. The output can be 
obtained either in frequency table or histogram form, the intervals for 
which are specified by the programmer. Mean values, standard devia- 
tions, maxima, minima and number of observations in each frequency 
class are provided as automatic output for table arguments. All attributes 
of a modeled system can conceivably be tabulated in a frequency table 
or plotted in a histogram. The format of a frequency table cannot be 
altered. The programmer can, however, program the dimensions of a 
histogram. 

3 . Comparison 

Data collection and display in FORTRAN is flexible, ob- 
trusive and programmer time consuming. In contrast, data collection 
in GPSS is non-flexible , unobtrusive and rapidly programmed. FORTRAN 
is flexible in that all data can be displayed according to any desired 
format. The GPSS format is compact and informative but it cannot be 
altered by the programmer. 

Data collection and display is obtrusive in FORTRAN because 
data collection methods can interfere and obscure the simulation logic 
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of the system. GPSS data collection and display statements do not inter- 
fere with program logic and can normally be placed in a separate portion 
of the program to avoid obscuring the logic of statement execution. 

FORTRAN data collection and display techniques require con- 
siderable programmer time in engineering the desired display format. 
Additionally, all desired output from FORTRAN programs must be pro- 
grammed since no output is provided automatically. GPSS automatically 
provides a considerable amount of output and allows the programmer to 
collect other data at his discretion. The programmer specified output 
is rapidly programmed within the simulation since the output format 
r does not have to be considered. 
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V. CONCLUSIONS 



The key words in describing the usefulness of FORTRAN for simulation 
is flexibility and universal availability. FORTRAN afforded flexibility in 
all the areas considered. The analyst may choose the most efficient 
methods to reproduce stocastic phenomena. The methods used to represent 
the static and dynamic behavior of the system were limited only by the 
ingenuity and resourcefulness of the analyst. This flexibility, however, 
may be gained at a stiff price. The price paid for flexibility was paid 
in man hours spent in conceptualizing the problem, programming, monitor- 
ing and debugging the FORTRAN model. For the inexperienced analyst, 
flexibility may spell trouble. It is possible for the analyst to select a 
comparatively inefficient method, since many diversified methods are 
available to model a given system. FORTRAN, regardless of any dis- 
advantages, will however still be used extensively for simulation. The 
reason being is its universal availability. FORTRAN compilers are 
available for nearly all commercially manufactured general purpose 
computers . 

The principal advantage of GPSS is that the language is user 
oriented. The GPSS language assists the programmer in conceptualizing 
the problem, and the GPSS model required reduced programming, debug- 
ging and monitoring time compared to the FORTRAN model. The GPSS 
block diagram flow chart closely represented the simulated system, and 
the flow chart block mnemonics were quickly convertible to the GPSS 
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computer programming language. The above contributed significantly to 
reducing the time required to obtain a working simulation model. The 
primary disadvantage of GPSS is that output is restricted to bar graphs 
and frequency tables. It is conceivable that the GPSS output might have 
to be reformated by clerical personnel to conform to specific report require- 
ments . 

The problem formulation, programming and debugging time that was 
saved by GPSS, in contrast to FORTRAN, was of the order of three or 
four to one. This time differential is of such magnitude that the analyst 
should consider the use of GPSS even if the analyst has had no previous 
experience with the language. The advantages of GPSS over FORTRAN 
for simulation indicate that GPSS should be used whenever both GPSS 
and FORTRAN compilers are available. 

It should be noted, however, that GPSS is basically restricted to 
the simulation of queueing systems while FORTRAN, as a general purpose 
language, can be used to simulate any system that can be described by 
a sequence of mathematical expressions. 
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APPENDIX A 



Exogenous Data 

The values utilized for exogenous data were hypothetical but 
conceptually realistic. One expects the mean length and arrival rate 
of messages to decrease as the precedence of the message increases. 
The actual mean values are normally a function of the type of military 
units within the geographical area serviced by the terminal. The arrival 
rate of messages is normally a function of the time of day. The arrival 
rate of messages generally peaks sometime after the end of the normal 
work day. Exogenous data that was utilized for both the FORTRAN and 
GPSS models is represented in Figures 5 thru 12 of this appendix. 
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PERCENT PEAK LOAD ARRIVAL RATE FOR TERMINAL ONE 
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PERCENT PEAK LOAD ARRIVAL RATE FOR TERMINAL TWO 
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Figure 6 



PERCENT PEAK LOAD ARRIVAL RATE FOR TERMINAL THREE 
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Figure 7 



PERCENT PEAK LOAD ARRIVAL RATE FOR TERMINAL FOUR 
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Figure 



PERCENT PEAK LOAD ARRIVAL RATE FOR TERMINAL FIVE 
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PERCENT TRAFFIC TRANSMITTED TO EACH TERMINAL 
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Figure 10 



MEAN PEAK HOUR ARRIVAL RATE BY PRECEDENCE 
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Figure 11 



MEAN MESSAGE TRANSMISSION TIME BY PRECEDENCE 
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Figure 12 



APPENDIX B 



FORTRAN Computer Model 
A. Definition of Arrays 

1. ARTE(4,5) contained the mean peak hour arrival rates listed 
in Figure 11, Appendix A. Columns corresponded to terminals 1 thru 5, 
respectively. 

2. BACKLP(IO) stored the maximum low precedence backlog of 
messages in minutes for each terminal and relay transmitter bay. 

3. BACKHP(IO) stored the maximum high precedence backlog of 
messages in minutes for each terminal and relay transmitter bay. 

4. FF(5,10,3) stored flash message parameters for the terminals 
and relay transmit bays. Rows 1 thru 5 stored arrival times, and position 
FF (1,1,1) contained the lowest arrival time for terminal one. Columns 

1 thru 10 represented terminal stations 1 thru 5, and relay transmitter 
bays 6 thru 10, respectively. Relay transmitter bay 6 transmitted messages 
to terminal 1, bay 7 transmitted messages to terminal 2, etc. The third 
dimension in the array was used to store message attributes as follows: 
Position 1 was arrival time, position 2 was message transmission time 
and position 3 was destination. 

5. HOLDl(i) - HOLD 10 (i) represented the circuits between 
terminals and relay transmitter bays. The mnmonics of the HOLD arrays 
corresponded to the numbering depicted in Figure 1, Section B, Chapter 
III. The dimension of each array was equal to the number of transmit 
circuits between two stations and was defined by variables HI thru H10. 
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6 . 



NEXEV(20,2) was the array used to store the next event times 



for the terminals, relay transmitter bays and groups of transmitter circuits. 
The transmitter circuits were represented by arrays HOLD1 thru HOLD 10. 

A message that was transmitted was stored in a HOLD array for the trans- 
mission length of the message. Column 1 of array NEXEV recorded the 
next event time and column 2 was a statement number corresponding to 
the subroutine call statement associated with that event. 

7. 00(20,10,3) stored immediate message parameters for the 
terminals and relay transmit bays. Rows 1 thru 20 stored arrival times 
and position 00(1,1,1) contained the lowest arrival time for terminal 
one. Columns 1 thru 10 represented terminal stations 1 thru 5, and 
relay transmitter bays 6 thru 10, respectively. Relay transmitter bay 6 
transmitted messages to terminal 1, bay 7 transmitted messages to 
terminal 2, etc. The third dimension in the array was used to store 
message attributes as follows. Position 1 was arrival time, position 2 
was message transmission time, and position 3 was destination. 

8. PP(30,10,3) stored priority messages for the terminals and 
relays. Rows 1 thru 30 stored arrival times and position PP (1,1,1) 
contained the lowest arrival time for terminal one. Columns 1 thru 10 
represented terminal stations 1 thru 5, and relay transmitter bays 6 thru 
10, respectively. Relay transmitter bay 6 transmitted messages to 
terminal 1, bay 7 transmitted messages to terminal 2, etc. The third 
dimension in the array was used to store message attributes as follows: 
Position 1 was arrival time, position 2 was message transmission time, 
and position 3 was destination. 
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9. RR (60,10,3) stored routine messages for the terminals and 
relays. Rows 1 thru 60 stored arrival times and position RR(1,1,1) con- 
tained the lowest arrival time for terminal one. Columns 1 thru 10 
represented terminal stations 1 thru 5, and relay transmitter bays 6 thru 
10, respectively. Relay transmitter bay 6 transmitted messages to 
terminal 1, bay 7 transmitted messages to terminal 2, etc. The third 
dimension in the array was used to store message attributes as follows: 
Position 1 was arrival time, position 2 was message transmission time, 
and position 3 was destination. 

10. ZLY(8 , 5) was initialized to contain the percent peak load 
for each station as represented in Figures 5 thru 9, Appendix A. 

11. ZLX(8 , 5) was initialized to contain the time of day cor- 
responding to percent peak load for stations 1 thru 5 as represented in 

T — 

Figures 5 thru 9, Appendix A. 

12. ZPTT(5 , 5) was initialized to contain the percent of messages 
transmitted among stations, as depicted by Figure 10, Appendix A. 

B . Subroutines 

1. BACKLOG (K) 

a. Arguments 

K was a terminal site number or relay transmitter bay 
number, as depicted by Figure 1, Section B, Chapter III. 

b. Purpose 

Subroutine BACKLOG maintained statistics on the maxi- 
mum high and low precedence backlog for each terminal and for each relay 



transmitter bay . 
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2. DEST(J, DES) 



a . Arguments 

J was a value 1,2 ,3 ,4 or 5 that represented a terminal 
designation. DES was a number 1 thru 5 corresponding to the destination 
terminal of the message. 

b. Purpose 

Subroutine DEST obtained a pseudo-random number and 
in conjunction with array ZPTT determined the destination of an originated 
message . 

3. FILL(K) 

a . Argument 

/ 

Argumertt K corresponded to the terminal designators 

1 , 2 , 3 , 4 and 5 . 

b. Purpose 

Subroutine FILL created a specified number of messages 
of each precedence to initially fill arrays RR, PP , OO and FF . Message 
transmission time and destination were also generated for each message. 

4. IHOLD(HOLD,L, M) 
a. Arguments 

HOLD was an array of dimension determined by variables 
HI thru H10. The dimension equaled the number of circuits between two 
stations. L was the dimension of array HOLD in the subroutine. M was 
the row number in array NEXEV(20,2) corresponding to the event under 
consideration . 
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b. Purpose 

Subroutine IHOLD was used to release a circuit to 
availability after a message was completely transmitted. If two events 
had equal next event times , subroutine IHOLD was called prior to sub- 
routines TERM or RELAY. 

5. INTER (K, TIME, ZNS) 

a. Arguments 

K was a value 1 thru 5 corresponding to a terminal 
designator. TIME was the time of day in seconds. ZNS was the percent 
of peak load corresponding to the argument TIME. 

b. Purpose 

Subroutine INTER performed a linear interpolation on 
Figures 1 thru 5 of Appendix A. INTER produced a percent peak load 
value for a given time of day. 

6. RANDO(IY, YFL) 

a. Arguments 

EX was in FORTRAN common with the main program and 
was initialized an odd integer number. After the initial value had been 
specified, EX became a function of IY which also was an integer. YFL 
was a pseudo-random number in the range (0,1). 

b. Purpose 

Subroutine RANDO generated pseudo-random numbers 

in the range (0,1). 
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7. 



RELAY (HOLD , K , L , M) 



a. Arguments 

HOLD was an array of dimension L. L was equal to 
the number of transmit circuits. K was the row number in array NEXEV 
(20,2) corresponding to the event under consideration. 

b. Purpose 

RELAY processed messages being transmitted from a 
relay transmitter bay to a terminal. The subroutine considered messages 
of precedence flash, immediate, priority, and routine on a first in first 
out discipline by precedence. Flash messages pre-empted a lower 
precedence message if no circuit was available. The message selected 
was transferred to the respective HOLD array where it remained for its 
transmission time. The backlogged messages at the terminal then 
advanced one position in their respective array columns. Necessary 
statistics were compiled to determine backlogs and circuit utilization. 
The next event time for the given terminal was set in accordance with 
circuit availability. 

8. RFILL(K , MOVE , PREC) 
a . Arguments 

K corresponded to the terminal station numbers 
1 , 2 , 3 , 4 , 5 and relay transmitted bay numbers 6,7,8,9,10. If MOVE 
equaled 1, RFILL was called by a relay and if 2, RFILL was called by 
a terminal . 
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Purpose 



b . 

Subroutine RFILL had two functions. For both the relay 
and terminal, the subroutine advanced quequed messages one position 
in the message’s respective precedence array column. After the messages 
had been advanced, at least the last position in the array was vacent. 

For a terminal, a new message was generated and was placed in the last 
position of the precedence array column. 

9. RELREO(K, PREC) 

a. Arguments 

K referred to relay transmit bay positions, and assumed 
values 6, 7, 8, 9 and 10. PREC assumed values of 1,2,3 and 4 that 
corresponded to precedences routine, priority, immediate and flash, 
respectively. 

b. Purpose 

Subroutine RELREO ordered message arrival times in 
relay transmit bay columns of the precedence arrays. This was done in 
order to assure that a message with the lowest arrival time was in 
position one for each transmit bay. 

10. TERM (HOLD, K,L,M) 
a . Arguments 

HOLD was an array of dimension L. L was the number 
of transmit circuits. K was the row number in array NEXEV (20,2) cor- 
responding to the event under consideration. 
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Purpose 



b . 

Subroutine TERM processed messages that were trans- 
mitted from a terminal to a relay. The subroutine considered messages 
of precedence flash, immediate, priority and routine on a first in first 
out discipline by precedence. Flash messages pre-empt a lower prece- 
dence message if no free circuit was available. A message was trans- 
ferred to its respective HOLD array and to a future events list in queue 
in the relay transmitter bay that served the desired destination. The 
backlogged messages at the terminal advanced one position, and then 
a new message was added to the future event list. Necessary statistics 
were compiled to determine backlogs and circuit utilization. The next 
event time for the given terminal was set in accordance with circuit 
availability . 

11. UTIL (CLOCK) 

a . Arguments 

Clock was a constant equal to the length of computer 
run in simulation time units, divided by 1,000. 

b. Purpose 

UTIL computed the average circuit utilization for each 
terminal and relay bay to three figures. 
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FORTRAN COMPUTER PROGRAM 
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DO 301 1=1, HI 
301 H0LD1 ( I )=90C0000 
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APPENDIX C 



GPSS Computer Model 



Appendix C contains a GPSS flow chart of the torn tape relay 



network and the GPSS computer program used to 



simulate this system. 
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GPSS FLOW CHART SYMBOLS 



(gen) 

BLOCK 

© 

<0 



GENERAT e 
TRANS ACTIONS 



4 3 GPSS BLOCly Types 
C ReP. 6 > 7 > 8 3 I 4 ) 



X I S A M M £ M 0 M I C. 
CONNECTOR 



OFF PAGE 
C OMMFCTOR 
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Gea yeR^re A Pouttaje 
MessAGi: for szre 



STORE RRTF 

OF PP*?I UAL FOR 
I\JS AT AAlrSS^Oe 

ASSIGAJ L£ AJ&TH OF 
G E FO P/I^PaO-. 

ETER 1 



TRQ wsf? r ro cire 



GdrJERATE A 

pRlORLry /^esSAG 6 

for srre 

Store mea/v r&te 
op arr iwaz. FtfR 
Neat message 

/4SS/&/0 LBNQrTH Op 

hessA oe to 

PARAA AFrTER. / 

TRA/os f izR To cite 
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GrlzAjeRATe A/ 

Xmms brR re 

Me^AG-e for sire 

Stdrf MfflA; aatf 

OP ARRIVAL TOR 
l\JEKT MESSAGE 



LEVOTH 
OF MESSAO-e TO 

Parameter I 



TRa^sfeR ro c/re 



Generate /A 
Plash messag-e 
for sire 

STOR^z Me/A/v RAre 

OF /^RRluAi. /TOR 

/vext Message 
ASS/G/v LE/vg-th 

OF MESS A 6 -E 7 ^ 
PAR am ETP R I 



trqvspbr ro cire 
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ASSIGN 



assign 



TRANSFER 




ASSIGN 



A SS IGN 



ASSIGN 



ASSIGN 



l 



ASSIGN b&STlNAT/ON 
TO PARAMETER 3 

ASS IGI\I pequ/Rek 
ROUTING- TO 
PARhNllBTB R H- 

TRQhs fbr mfs' S. 
PRO/a C/7ES /j 2^3 Vj 5* 

TO xv\ir 



A S 5 IGN Rour/AOr 
FOR Z 57 HOP TO 
PARAMETER 5 " 

ASS/ 6 - A/ LO\ajBST 

XMTR. A/O, FOR 

TA/s sire ro 

Parameter 8 

ASS/G-AZ HIGHEST 
XfATR A/O FOR 
This s/re to 

P-aKfl K P TF J ^— 9 

pp PR EC. - 2 .J 00 - 3 

pP ; Vj RR~ JT due/ / 

RESER^b FOR P/?eE*p 7 Fr> 



128 



A 




ASSIGN USER CHAIM 

Nu^QbR "For bach 
5 /re to Par n iz 

SetoO HIGH PREC. 
he^SPOES TO HIRRI 



/QCCcmuL47"E LO\jJ 

Precedence mes. 
Backlog- statistics 



A SAveu/uwe 
LOCATION TO 
PAM. IO 

IP current* lovj PRsc. 
BRCPLOG* Less THAU 
PReu/Owi MA^ 
TRAuSPBR TO SLCTI 

Record Current 
Backlog- as 

MAhKIKu^ 

rRANSFfiR TO 
Slot / 
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SAi/ev/flLue 



ASS/G-M 



ASS/OM 



TEST 



S^vevALMe 



TRANSFER 



ACCUmuLA^ WIG-H 
PReciz DBfijce Ales, 
Backlog- srAr/sr/cs 



A 5Al/i?i/ALC^c 
Locat/om TO 
pflRV\. IO 



A SA^ei/ALue 

LOCAT/OAj to 

PARM. £ 



IP Cu RRCjoTHIGW PR&C. 
BACKLOG- Less TRAM 
Previous max 
TRAMS F6R TO StCTj 



ReCORto Cui*R(rMT 
&ACKLO G- AS 
M AX I lAUM 



TRAM S F (r R TO 

Slc ri 
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1CTI 




Se^ecr A i^Ree XhrR. 
IF move Available 
TRANSFER 7b flash 

TRAK)^>fer to 
Available am tR. 

Go ro trakj^ 



^ESSACrF /s stored 

JM USER CHRI/J UWTH 
TIME ROR XK1^S7 0A) 



IF hvrSSAGe A/OT 
R-As/^G-oro user 
CH ATM ( LIKJHI ) 

select the x mt r< 

■S0AJ Dl NJ (> LOCuEST 
PR0C fr DEM C E MBSS, 

^LASH Mess AGE PW-Bntr 
Loweyr PRbo 
PRB-Bh^T^O mbssagg 
G-O TO U SB R CHR7/0 

teCRBiAeAJT H/6/4 
PR6C. &AC/O-0G 13V 
TTMB LENGTH 
MBSS AG fr 
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A 



flDV/AwClz 



return! 



Uajunk 



TRAft/SFeR 



Hold XmtR por 

L6A/G -TH OP 

wessRG-e 



Returm XmtR por 
Loujer Precede/vce 
MESSAO-E s 

UN LIN K Q MESS A^c- 
FROf^ USER QHAI/V 
TO BNTBR M 5LC rr 

tRRwsPe/? td resTL 
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5 ei 2 r£ A XmR.r/ye 

AjUl^BER PCCOR bl^Or 

TO 1//9LU0 op 





S/H/ei/Ai-ue 








(^!)«« — 


res r 











resr 










SAuei/ALue 








(wvj)*f — 


.transfer 



RecoR o Pppc 

OP A'UrSSAGIc HOLD* MG- 
ErACH XMTR. 

IP Messooir WPS 
PRlrVIOUSV PRe-EnPTb 
TRA/o^FeR 



IF h 0 SSAOIr /S 
PRIORITY OR ROUT/AJE 
IROAJ SF F R O P/U ” 



Dec rpkfmt h/^h 

PR EC 0 ?Df ^CF /3AC/U06 



TRfWSFlSR TO A&VI 
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UDPSI 




DecRerAEMT lo PRec 
B^ctaoc- 8V niAE 

leMGTH FLAS H M£SS. 



AOV/A/^CG 



Hold Xmtr For 

LG/O&TH OF Aies.s. 





Reie, 


as £ 


amt/R. 

AjUMQCR Z)ccorMiu6- 

TO 7 










UN LIN k 


UNLINK 4 MESS40-B 
PRON\ U^&R CHQIV 
TO B/oTBR BJ-OC/C 






" %cr/ u 


>- — 


trapper, 


TRA ajSFCR To 
’T^ST z " 
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XF WESSAO- E OV 2 ncJ 
Le 6 - OF ROUTJjoG- 
TRA/oSFeR TO LECrZ 



/ 3 SS/ 6 AJ ROUT^Cr- 

pOR 2.0^ /VoA 



>^T COOWIF/C 
pc AJUfAtifr^ CP 



XF AAESSAGf OAJ 
|SJ /Vo PjTRA/USFeR 

To lopf , 4. e. Rspbat 

eve /e 

7 "/r PNM ATI: 

O AJ 2 ^ NO & 
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GrevtRQTZ T KANSACT/oAJ 
LjITH 3 0 PfiRMZT5R% 




AWAKldz 



AS SfG /O 



ft SSI 6 At 



ASSIGN 



TQ&ULRTft 



TA BLf= 



TE - Ry \ i *> p)TF 

i 



iJOLCs TR A A; ^ACTJoaj 
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